Wavepacket approach to the cumulative reaction probability within the flux operator formalism.

Expressions for the singular flux operator eigenfunctions and eigenvalues are given in terms of the Dirac delta-function representable as a localized Gaussian wavepacket. This functional form enables computation of the cumulative reaction probability N(E) from the wavepacket time-correlation functions. The Gaussian based form of the flux eigenfunctions, which is not tied to a finite basis of a quantum-mechanical calculation, is particularly useful for approximate calculation of N(E) with the trajectory based wavepacket propagation techniques. Numerical illustration is given for the Eckart barrier using the conventional quantum-mechanical propagation and the quantum trajectory dynamics with the approximate quantum potential. N(E) converges with respect to the Gaussian width parameter, and the convergence is faster at low energy. The approximate trajectory calculation overestimates tunneling in the low energy regime, but gives a significant improvement over the parabolic estimate of the tunneling probability.